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SUBSONIC,  TRAN'S  ONTO , 


SUPERSONIC  NOZZLE  FLOW  FIELDS 


Lin  Tung  Chi  and  Chia  Chen  Hr.ueh 
ABSTRACT 

The  present  article  gives  the  general  form  of  the  equation 
for  the  two-dimensional  constant  steady  irrotational  isentropic 
flow  of  an  ideal  gas  in  conformal  curvilinear  coordinates.  Sim¬ 
plified  equation  for  the  stream  function  and  its  general  solu¬ 
tion  are  given  using  the  streamlines  and  equipotential  lines  of 
the  corresponding  incompressible  potential  flow  as  coordinates. 

The  results  are  applied  to  nozzle  flow,  giving  the  solution  for 
nozzle  flow  from  subsonic  through  transonic  to  the  supersonic 
range,  for  any  desired  choice  of  radius  of  curvature  of  throat 
wall,  contraction  ratio  and  largest  angle  of  inclination  of  the 
wall.  This  solution  is  valid  for  different  ratios  of  specific 
heats . 

As  an  example  of  the  application,  we  have  calculated  the 
flow  characteristics  of  a  typical  nozzle.  These  include:  line 
of  constant  Mach  number  for  subsonic,  transonic  and  supersonic 
nozzle  flow;  velocity  line,  line  of  constant  Mach  number,  influence 
line,  limiting  characteristic  line,  branchline  and  equit.emporal 
lines  for  supersonic  nozzle  flow. 

This  method  can  be  extended  to  flow  around  an  obstacle,  over- 
weir  flow  in  the  nozzle  and  flow  through  a  grating.  Particularly, 
one  can  obtain  better  results  In  the  transonic  region.  Further¬ 
more,  the  method  can  be  extended  to  cases  where  there  are  chemi¬ 
cal  reactions  under  equilibrium  or  non-equilibrium  conditions, 
or  where  there  exists  axial  symmetry. 


SYMBOLS 


function.  Equation  (14) 

a  velocity  of  sound, 

Equation  (3) 

B,  (£,a)  function.  Equations 
(33),  (3*0 

b  half-height  of  nozzle, 

throat.  Equation  (1) 

D  determinant.  Equation  (9) 

Er.  (  £  ,n  )  function  relating  to 
the  transformations. 
Equation  (14) 

f(C)  transformation  of  confor¬ 
mal  curvilinear  coordi¬ 
nates,  Equation  (6) 

G1(£,n)  function.  Equation  (.33) 

g(£)  flow-related  function. 
Equation  (12) 

gg(0)  flow-related  function  for 
maximum  flow  at  the 
throat.  Equations  (18), 

(.22  ) 

H  (£,n)  amplification  factor 

in  going  from  z-plane  to 
C-plane,  Equation  (6) 

1  point  of  inflection  at  the  wall 
Figure  1 

i  =  /^T 

k  =  nozzle  flow  factor. 

Equation  (25) 

M  Mach  number.  Equation  (4) 

N  (£,n,ct)  functions. 

Equations  (33),  (34) 

n  nozzle  contraction  ratio 
Equations  (31),  (32) 

p  pressure,  Equation  (3) 

q  velocity.  Equation  (1) 

r  radius  of  curvature  of  throat 
wall,  Equations  (31),  (32) 

s  distance  along  flow  line, 
Equation  (27) 


T  absolute  temperature  of  gas, 
Equation  (2) 

u  velocity  component  along  £ 
direction.  Equation  (2) 

v  velocity  component  along  n 
direction,  Equation  (2) 

x  rectangular  coordinate. 

Equation  (6) 

y  rectangular  coordinate, 

Equation  (6) 

z  =  x  +  iy  complex  variable  in 
rectangular  coordinates. 
Equation  (6) 

a  =  ch-^  (ochX),  Equation  (33) 

8  =-iot,  Equation  (37) 

Y  ratio  of  specific  beats.  Equa¬ 
tion  (3) 

6  =-tg-1  )  angle  between 

flow  and  n  direction. 
Equation  (26) 

C  =  £  +in  complex  variable  in 

conformal  curvilinear  coor¬ 
dinates,  Equation  (6) 

q  conformal  curvilinear  coordi¬ 
nate,  Equation  (6) 

rih  nozzle  shape  factor. 

Equation  (28) 

9,  largest  angle  of  inclination 
of  nozzle  wall.  Equations 
(3D,  (32) 

X  nozzle  shape  factor.  Equation 
(.28) 

y  =  tg-1  (M2-l )-1//2  Mach  angle. 
Equation  (26) 

v  positive  whole  number.  Equation 

(14) 

£  conformal  curvilinear  coordi¬ 
nate,  (Equation  (6) 

£,  point  of  inflection  at  wall, 
Equations  (31),  (32) 
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p  density  of  pas,  Equation  (3) 
o  nozzle  shape  factor.  Equation  (28) 

t  time  spent  in  nozzle  by  ras  particles,  Equation  (27) 

<t>  potential  function.  Equation  (5) 

ip  stream  function,  Equation  (5) 

i(i  largest  flow  in  nozzle.  Equations  (19),  (?3) 

S 

w(£)  function.  Equation  (25) 

Subscripts : 

*  critical  parameter  for  which  stream  velocity  equals  local 
velocity  of  sound 

b  nozzle  wall  parameter 
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ordinal 
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partial 
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INTRODUCTION 


Nozzles  have  wide  applications  in  engineering.  As  early  as 
the  80's  of  the  19th  Century,  Laval  [1]  used  convergent-divergent 
nozzles  in  steam  turbines  to  obtain  supersonic  flows.  During 
the  same  period  of  time,  Reynolds  [2]  analyzed  theoreitcally  one- 
dimensional  supersonic  nozzle  flows.  Early  in  this  century  Meyer 
[3]  was  the  first  to  analyze  two-dimensional  supersonic  nozzle 
flows.  In  the  30's,  Taylor  [4],  Hocker  [5]  and  Goertler  [6] 
separately  studied  the  higher  subsonic  nozzle  flow  near  the  thorat 
wall  that  is  partially  supersonic.  After  the  40's,  the  develop¬ 
ment  of  supersonic  aircraft  and  jet  propulsion  technology  spurred 
research  on  supersonic  flows,  resulting  in  a  large  quantity  of 
work  on  theoretical  analysis  and  designing.  From  the  viewpoint 
of  applications,  when  the  ratio  of  the  radius  of  curvature  of  throat 
wall  to  the  half-height  of  the  throat  passage  is  very  large,  the 
design  problem  of  the  supersonic  nozzle  is  basically  solved.  In 
the  past  decade,  high  acceleration  nozzles  became  necessary  to 
shorten  the  nozzle  length  of  the  generator  of  rockets  and  to  in¬ 
crease  the  efficiency  with  which  the  gas-propelled,  excited  light 
nozzle  freezes  high  temperature  gases  under  non-equilibrium  condi¬ 
tions.  In  order  to  achieve  this,  the  radius  of  curvature  of  the 
throat  wall  needs  to  be  decreased  and  the  angles  of  convergence 
and  divergence  have  to  be  increased.  The  presently  available  anal¬ 
yses  and  design  methods  for  supersonic  nozzle  flows  do  not  meet 
this  design  requirement  and  hence,  are  not  applicable  to  the  study 
of  supersonic  flows  in  nozzles  with  small  radius  of  curvature  of 
throat  wall. 

A  large  quantity  of  information  has  been  accumulated  regard¬ 
ing  the  study  of  supersonic  nozzle  flows.  According  to  methods  of 
treating  the  problem,  one  can  divide  it  into  eight  classes: 
expansion  in  power  series  [3-6,8,15,51],  hodograph  method  [7,9,10, 
19,23,25,27,32734,37,44 ,46-48,50,52-58,60] ,  iterative  substitution 
relaxation  method  [13,14],  inversion  method  [16,43,67-69,80,87], 
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expansion  In  snail  parameters  [59,64,65,70,93],  band  interration 
[61,62,63],  time  correlation  method  [71,72,77,79,81,84,85]  and 
method  of  imaginary  characteristic  line  Prandtl's  relation  [82]. 


In  the  method  of  expansion  in  power  series,  the  velocity 
potential  or  velocity  component  is  expressed  as  a  power  series  in 
the  position  coordinates.  A  definite  velocity  distribution  is 
assumed  and  an  approximate  solution  for  the  nozzle  flow  is  ob¬ 
tained  by  finding  the  expansion  coefficients.  For  example, 
Oswatitsch  et  al»[8]  assume  that  all  the  equivelocity  lines  are 
quadratic  curves  in  positive  coordinates.  Sauer  [15]  assumes  that 
the  velocity  distribution  on  the  central  streamline  at  the  throat 
bears  a  linear  relation  to  the  abscissa.  This  method  gives  an 
approximate  solution  for  the  case  of  larger  radius  of  curvature 
of  throat  wall. 

In  the  hodograph  method,  the  velocity  components  are  used  as 
coordinates,  and  the  nozzle  flow  properties  are  studied  on  the 
hodograph  plane.  Accordingly,  the  equations  become  linear,  but 
the  boundary  conditions  become  more  complicated  [91].  There  are 
many  papers  in  which  the  nozzle  flow  is  treated  using  the  hodograph 
method.  Ringleb  [7]  was  the  first  to  use  this  method  to  obtain 
the  singular  solution  -for  flow  around  a  semi-infinite  plane  wall 
that  contains  partial  supersonic  regions.  Temple  et  al.  [9]  eluci¬ 
dated  this  flow  and  pointed  out  that  it  is  equivalent  to  the  nozzle 
flow  near  the  throat  wall  that  is  partially  supersonic.  Frankl  [19J 
and  Lighthill  [27]  independently  introduced  the  concent  of  branch¬ 
line  of  transonic  nozzle  flow  in  the  hodograph  plane  and  used 
methods  of  approximation  to  discuss  the  flow  characterics  in  the 
vicinity  of  the  velocity  point  on  the  center  line.  Tomotika  et  al . 
[32]  used  an  assumed  gas  to  analyze  changes  in  subsonic  and  super¬ 
sonic  nozzle  flow.  This  has  definite  significance  in  the  under¬ 
standing  of  flow  characteristics  in  the  throat  region.  Cherry 
[34,37,52,56]  used  the  method  of  iterative  addition  of  potentials 
to  obtain  an  accurate  solution  for  the  supersonic  flow  in  nozzles 
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with  particular  wall  contours.  This  study  is  important  in  con¬ 
tributing  to  the  understanding  of  the  flow  in  the  transonic 
region  of  the  nozzle.  Because  the  boundary  conditions  are  diff¬ 
icult  to  control,  it  cannot  be  conveniently  applied  in  engineer¬ 
ing  designs. 

The  iterative  substitution  relaxation  method  is  based  on  the 
relaxation  method.  The  supersonic  nozzle  flow  is  discussed  using 
incompressible  potential  flow  and  eauipotential  line  as  coordi¬ 
nates.  Green  et  al.  [13]  gave  the  basic  equations  for  this  method 
and  Pox  et  al- [14]  applied  the  method  to  an  example  of  supersonic 
nozzle  flow  problems.  Because  of  the  limitations  on  the  choice 
of  the  form  of  the  solution  and  the  boundary  conditions,  no  prac¬ 
tically  useful  results  were  obtained. 

The  inversion  method  was  first  developed  by  Friedrichs  [lo] 
who  also  treated  the  problem  of  nozzle  flow.  Liepman  [43.]  used 
this  method  to  treat  two-dimensional  nozzle  flow;  Hopkins  et  al. 
[67,69]  and  Pirumov  et  al  [68]  separately  treated  axially  symme¬ 
trical  nozzle  flow.  This  method  converges  rather  slowly  and  the 
boundary  conditions  cannot  be  effectively  controlled.  Brodsky 
[78]  adopted  the  method  of  finite  differences  to  solve  the  equa¬ 
tions  of  the  inversion  method  and  discussed  supersonic  nozzle  flow. 
He  found  that  oscillations  of  the  solution  occurred  in  the  sub¬ 
sonic  region.  The  reason  was  that  the  boundary  conditions  for 
elliptical  equations  could  not  be  determined  from  the  velocity 
distribution  on  the  center  line  alone. 

Expansion  in  small  parameters  was  first  developed  by  Hall  [59] 
who  used  the  reciprocal  of  the  radius  of  curvature  of  the  throat 
wall  as  the  small  parameter  and  expressed  the  geometric  boundary 
of  the  nozzle  and  the  velocity  components  in  terms  of  power  series 
in  this  small  parameter.  Assuming  a  given  relation  between  velo¬ 
city  and  the  Abscissa,  he  obtained  the  first  three  orders  of  the 
solution.  The  first  order  solution  was  obtained  under  the 
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assumption  o  1'  a  linear  relationship  between  the  velocity  and  the 

abscissa.  The  result  arrees  with  Sauer's  [15].  The  second  and 

third  order  solutions  are  improved  over  those  obtained  by  Sauer 

~~  -1  1  1  !  .  ■ ) 

[15].  Similarly,  Kliegel  et  al.  L70]Vusod  the  r-cTprocal  of  the 
radius  of  curvature  of  the  throat  wall  as  the  small  parameter.  He 
found  an  approximate  solution  to  the  entire  nozzle  flow  by  includ¬ 
ing  in  the  one-dimensional  flow  a  correction  in  the  small  para¬ 
meter.  The  above  results  are  applicable  to  the  case  of  large 
radius  of  curvature  of  throat  wall.  Kliegel  et  al.[73]  used 
double  cylindrical  coordinates  and  took  the  reciprocal  of  radius 
of  curvature  of  throat  wall  increased  by  1  as  the  small  parameter 
in  the  expansion  of  the  velocity  components.  By  comparing  with 
Hall's  [59]  series,  the  coefficients  of  expansion  were  determined. 
This  improved  the  convergence  of  Hall's  series. 

Band  integration  is  a  method  that  employs  band  integration  of 
supersonic  flow  around  a  blunt-headed  body  to  calculate  nozzle 
flow.  Alikhashchkin  et  al.  [61]  used  one  and  two  bands  separately 
to  compute  the  flow  characteristics  of  a  two-dimensional  nozzle 
with  given  shapes  of  the  wall.  Holt  [62]  used  one  band  to  compute 
the  surface  contour  lines  of  a  nozzle  with  given  velocity  distri¬ 
bution  on  the  center  line.  Favorskiy  [63]  used  two  bands  to  com¬ 
pute  the  flow  characteristics  of  an  axially  symmetric  nozzle  with 
given  wall  contours.  In  this  method,  each  band  contains  a  "saddle 
point"  as  a  singular  point  which  places  a  limit  on  the  available 
number  of  bands  in  practical  applications.  Besides,  the  method 
integration  only  gives  the  average  value  between  the  bands  and 
cannot  be  used  to  characterize  local  flow.  In  general,  this  method 
is  not  suitable  for  nozzles  with  small  radius  of  curvature. 

In  the  time  correlation  method  for  the  solution  of  transonic 
flow,  the  mixed  equations  are  reorganized  into  hyperbolic  equa¬ 
tions  by  introducing  the  time  variable.  The  method  of  differences 
is  used  to  obtain  a  solution  for  unstable  motion  of  the  gas  des¬ 
cribed  by  the  hyperbolic  equations.  In  the  limit  of  infinite  time, 
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an  approximate  solution  for  stable  flow  is  obtained.  This  method 
offers  a  ne.-:  approach  to  the*  solution  of  the  transonic  flow  prob¬ 
lem.  This  method  is  one  of  the  most  effective  for  non-isent ropic 
flows,  i .  e . ,  when  there  are  excited  waves.  The  major  works  using 
this  method  are  by  Migdal  et  al.  [71],  Ivanov  et  al.[72],  Wehofer 
et  al. [77],  Laval  [7*0,  Serra  [8l],  Brown  et  al. [84],  Cline  et  al. 
[85].  Because  the  time  variable  is  introduced  to  convert  a  con¬ 
stant,  stable  flow  problem  into  one  with  unstable  flow,  the  prob¬ 
lem  is  complicated  by  the  additional  independent  variable. 


The  method  of  imaginary  characteristic  line  Prandtl's  rela¬ 
tion  came  about  as  follows.  Armitage  [66]  and  Rao  et  al. [7*0  used 
the  method  of  imaginary  characteristic  line  to  compute  the  flow 
in  the  subsonic  section.  Extension  to  the  sound  velocity  line 
was  made  and  Prandtl’s  relation  was  used  to  treat  the  so-called 
mixed  region  between  the  sound  velocity  line  and  the  limiting  char¬ 
acteristic  line  [82].  Writing  the  series  solution  for  the  mixed 
region  in  the  form  of  Prandtl's  relation  brings  about  convergence 
only  in  the  form  of  the  solution,  the  actual  solution  of  the  sub¬ 
sonic  flow  region  remains  complicated.  Furthermore,  there  is  also 
the  problem  of  matching  the  solutions  in  the  two  regions. 

In  summary,  there  are  many  methods  for  studying  nozzle  flow. 
However,  the  properties  of  the  mathematical  equations  for  the 
transonic  flow  region  at  the  throat  presented  great  difficulties 
in  the  study  of  nozzle  flow  and  some  of  the  important  problems  of 
nozzle  flow  do  not  have  a  satisfactory  solution.  For  instance, 
the  relationship  between  the  flow  characteristics  in  the  transonic 
region  at  the  throat  and  the  shape  of  the  wall  are  difficult  to 
determine,  especially  for  very  small  radius  of  curvature  of  the 
throat,  in  which  case  the  flow  in  the  throat  region  is  very  non- 
uniform.  Moreover,  the  lack  of  a  simple  method  to  solve  the  prob¬ 
lem  of  transonic  flow  makes  it  very  difficult  to  determine  the 
initial  conditions  in  the  design  of  the  divergent  section  of  the 
nozzle . 
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Owing  to  the  above  reasons,  empirical  or  semi-empirical 
methods  have  to  be  used  in  the  engineer  in."  dvr.  I  ,-,n  of  nozzles  at 
present.  Those  include  the  follow! nr,:  The  method  of  Busemann- 
Atkin  [16,28,3.1,88,92]  converts  radial  flow  into  parallel  flow. 
Because  of  its  dependence  on  design  experience  and  skills  in  the 
determination  of  throat  wall  contour,  this  method  has  consider¬ 
able  variabi 1 ity .  Puckett  et  al.  [ 24 , 29 , 30 , 83 , 86 ]  used  the  assump¬ 
tion  of  linear  velocity  line  in  specifying  the  initial  conditions 
for  the  design  of  the  divergent  section.  This  assumption  applies 
approximately  in  the  case  of  large  radius  of  curvature  of  throat. 

For  the  case  of  small  radius  of  curvature  of  throat,  where  the  velo¬ 
city  line  is  far  from  being' linear  [75,76]  this  assumption  is  not 
reasonable  [82,83].  Others  have  used  the  approximate  results  ob¬ 
tained  from  the  above  theoretical  analyses  as  initial  conditions  in 
the  design  of  the  divergent  section.  Examples  are  found  in  the 
series  method  of  Sauer  [15],  expansion  in  small  parameters  by  Hall 
[59],  and  other  semi-empirical  methods  [38,41].  In  the  design  of 
the  supersonic  section,  Guderly-Rao  et  al. [40 ,42 , 4 5, 49]  designed 
large-thrust  nozzles  with  given  initial  conditions  according  to  the 
principles  of  component  transformation.  At  present,  the  design 
methods  used  in  engineering  apply  only  to  nozzles  with  large  radius 
of  curvature  of  the  throat  wall.  The  problem  of  initial  conditions 
for  the  case  of  small  radius  of  curvature  is  what  needs  to  be 
solved  urgently. 

With  respect  to  the  above  problem,  we  have  given  the  general 
form  of  the  equation  for  the  two-dimensional  constant  steady  irro- 
tational  isontropic  flow  of  an  ideal  gas  in  conformal  curvilinear 
coordinates.  Simplified  equation  for  the  stream  function  and  its 
general  solution  are  given  using  the  streamlines  and  equipotential 
lines  of  the  corresponding  incompressible  potential  flow  as  coor¬ 
dinates.  The  results  are  applied  to  nozzle  flow,  giving  the  solu¬ 
tion  for  nozzle  flow  from  subsonic  through  transonic  to  the  super¬ 
sonic  rang/'  with  any  desired  choice  of  radius  of  curvature  of  throat 
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wall,  contraction  ratio  and  largest  ancle  of  Inclination  of  the 
wall.  As  an  example,  we  have  calculated  the  flow  characteristics 
of  a  typical  nozzle. 


II.  FLOW  EQUATION'S  IN  CONFORMAL  CURVILINEAR  COORDINATES 

Consider  the  two-dimensional  constant  steady  irrotational 
isentropic  flow  of  an  ideal  gas.  Dimensionless  quantities  are 
introduced  according  to  the  symbols  listed  in  the  list  of  symbols. 
These  are  denoted  with  a  above  the  corresponding  symbol: 


x/b. 

y  —  y/b. 

I  —  x  +  iy , 

r  —  r/b 

«/«*» 

v  —  v/a„ 

q  —  (u*  +  i»2)f 

»/»*, 

P  “  p/p*. 

A  —  p/p*,  — 

f  -  T/T, 

rajb. 

#  “  <P/a,b, 

mm  il,/a,p*bt 

M  ■**  q/a 

In  this  paper,  only  dimensionless  quantities  will  be  used. 
Hence,  the  will  be  omitted. 


The  conservation  of  mass,  momentum  and  energy  of  an  ideal  gas 
in  two-dimensional  constant  steady  irrotational  isentropic  flow 
can  be  represented  in  arbitrary  orthogonal  curvilinear  coordinates 
(Cjh),  respectively,  as: 


i  +  (Hipv)'  —  0 

—  (W,?)4  “  0 


(2) 


in  which  W,  —  -t-  y?)*  are  the  reciprocals  of 
Lamd  coefficients  in  the  orthogonal  curvilinear  coordinates  (£,n) 
[89]. 


The  state  equation,  isentropic  relation  and  the  relation  be¬ 
tween  sound  velocity  a  and  temperature  T  of  an  ideal  gas  are  given 

I 

respectively  by: 

P~pT,  a-  T*  (3) 
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From  Equations  (2)  arid  (3),  one  obtains  the  relation;-,  among 
density  p,  Mach  number  M  and  velocity  q  to  be: 


p 

M 


JL±±  _  JLHi.  9»YHt 
.  2  2  / 


O) 


From  the  equation  for  conservation  of  mass,  Equation  ( Pa  ) ^  ^ 
and  that  for  conservation  of  momentum,  i.e.,  irrotationalit.y , 
Equation  (2b),  we  define  stream  function  4>  and  potential  function 
<p  to  be  respectively: 


<J>t  *“  —pH\u,  </>,  “  pH,u 

C>t  •—  H,u,  «*  Htv 


(5) 


When  i.e.,  when  the  transformation  from  orthogonal 

linear  coordinates  x,y  to  curvilinear  coordinates  £,n  is  conformal, 
according  to  the  theory  of  complex  variable  [90],  the  following 
relation  exists  between  coordinates,  x,y  and  5,q: 


*-/<{>.  * 
H,  -  Ht  -  H 


*  +  iy ,  C 

ds 


5  +  *n 


(O 


where  H  is  the  amplification  factor  of  z  plane  with  respect  to 
[w.  r.  t.]  s  plane.  In  conformal  curvilinear  coordinates,  from 
equations  (4a)  and  (5),  one  obtains  the  following  relation  between 
density  p,  velocities  u,  v  and  stream  function  \p. 


pi-H  _  JL±1  pi  +  o 

u  *“  ,  •>  —  —<l>i/pH 


(7) 


Using  H 1  =  H 2  =  H  and  Equation  (4a),  one  can  eliminate  u,v,  <t> 

<J>^  and  p^,  from  Equation  (5)  to  obtain  the  equation  for  stream 
function  in  conformal  curvilinear  coordinates 


1) 


The  letter 
the  set. 


in  parentheses  denotes  the  order  of  an  equation  in 
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WCW-V  -  </>’)</>«  +  2W, +  H(WVpl  -  <A})<A„ 

-  (</-}  +  tf,<A,)  -  0  (8) 

The  determinant  for"  Equation  (8)  is 

D*  -  -4HVp*(l  -  A/1)  (a$ 

when  M  <  1,  Equation  (8)  becomes  an  elliptical  equation;  for  M  > 

1,  it  is  a  hyperbolic  equation. 

When  H=l,  Equation  (8)  can  be  simplified  to  obtain  the  equa¬ 
tion  for  the  stream  function  in  orthogonal  linear  coordinates  [9**] 

(«V  —  +  2  4>,4>i4>„  +  («V  —  <!>\  )4*„  —  o  (io> 

The  flow  equation.  Equation  (8),  given  in  this  paper  is  valid 
in  any  conformal  curvilinear  coordinate  system.  The  different  con¬ 
formal  curvilinear  coordinates  will  only  differ  in  the  expression 
for  the  amplification  factor  H  in  Equation  (8).  Thus,  the  compli¬ 
cated  relation  between  the  coordinates  is  included  into  the  ampli¬ 
fication  factor  H,  allowing  more  complicated  relations  between  the 
coordinates  to  be  used  in  practical  applications.  The  nozzle  coor¬ 
dinates  used  in  Section  V  are  an  example.  Compared  to  the  poten¬ 
tial  function,  the  advantage  of  using  the  equation  for  the  stream 
function  is  that  the  stream  function  is  a  constant  at  the  boundar¬ 
ies  and  is  thus  easier  to  treat.  We  will  discuss  this  in  greater 
detail  in  the  paper. 

III.  THE  SIMPLIFIED  FLOW  EQUATION  AND  ITS  SOLUTION 

In  two-dimensional  flows,  take  the  incompressible  flow  line 
and  the  corresponding  equipotent ial  line  as  curvilinear  coordi¬ 
nates.  In  such  coordinate  systems,  u>>v,  or  !  ^  !  >:>  I  ^  I  •  Accord¬ 
ing  to  Equation  (8),  after  elimination  of  and  Its  derivatives 
w.r.t.C.n,  one  obtains  the  simplified  flow  equation 

HVPv„  -  //,*;  -  o  (n) 

12 


r 


Equal  Ion  (11)  1 a  nonlinear  second  order  ordinary  di  fforential 

2 

equat  ion  in  n ,  wit.h  a  paramet  ric  variable  £ .  Eliminate  i !>r  from 

o  a 

Equation  (7a)  to  solve  for  H"  .  Substitute  this  into  Equation  (11) 
to  obtain  the  first  integral  of  the  simplified  flow  equation, 

( Equation  (11),  as 

<l>n  —  KOp(I>  t)  (12) 

in  which  g(  £ )  is  related  to  the  property  of  flow  and  the  amount  of 
flow.  It  is  a  function  of  £  and  is  called  flow-related  function 
for  simplicity.  In  flow  around  an  object,  it  is  related  to  the 
Mach  number  of  the  oncoming  stream;  in  internal  flows,  it  is  con¬ 
nected  to  the  amount  of  flow  passing  through  the  duct.  g(£)  is 
determined  from  boundary  conditions. 


Using  n  =  0  as  zero  streamline,  i.e.,  when  t|/(£, 0),  the  general 
solution  for  i|/(£,n)  can  be  obtained  from  Equations  (ba)  and  (12): 


t<{)  J;[i 


dn 

r  +  1  HJ(f,t;)J 


(13) 


For  a  given  conformal  curvilinear  coordinate  system,  the  ampli¬ 
fication  factor  H(£,n)  is  a  known  function  of  £,n.  In  Equation  (.13), 
only  the  flow-related  function  g(£)  is  an  unknown.  It  is  deter¬ 
mined  by  boundary  conditions  and  will  be  discussed  in  the  next  sec¬ 
tion.  The  solution  for  i|>(£,n)  given  in  Equation  (.13)  is  valid  for 
any  value  of  the  ratio  of  specific  heats  y 


Let  y  =  1  +  When  v  is  a  positive  whole  number,  substitute 
this  expression  into  Equation  (13)  and  expand  the  function  to  be 
integrated.  Equation  (13)  can  then  be  written  as  ( 


_  w 

*(5,»?)  -  2  '^E*(?.tVl'M(5),  E«(f,  .,)  -  ('  rj^dr) 

»**o  'it 

/*«-(-!)«  ci-±v(v-  l)...(v~*+  1) 

V-  1,  2,  3,  •••;  *-0,  1,  2,  v  >  \ 


(M) 
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For  a  given  set  of  conformal  curvilinear  coordinates,  £,  ,n, 

Ek  is  n  known  function  of  the  coordinates.  For  convenience, 

we  call  it  the  function  relating  the  transformations.  is  a 

function  of  y  only.  Thus,  the  stream  function  ']<(£, n)  becomes  a 
polynomial  in  g(£)  of  the  (2v  +  l)th  degree,  that  has  coefficients 
products  of  and  E^  (£,n)  and  that  contains  only  the  odd-powered 
terms . 

For  convenience  in  applications,  we  list  the  following  express¬ 
ions  for  t|/(£,n)  for  several  values  of  y. 

r  “.2  *(f,  fi)  -  Et(s,  v)g(S)  -  t 

Z  2 

y  -  1-5  <K|,  t?)-^  E.(f,  -  t  Et(S,  nV(?)  +  ~  * 

10  8  16 

ao  (15> 

r  -  1.33  *(g,  ,)  -  222  Sods,  rigdO  ~  JJ  £,(£,  ,)*»(*) 

216  72 

+  ^  EiCf ,  «?)*’  (£>  -  ^  £»(?» 


IV.  MAXIMUM  FLOW  AND  FLOW  CHARACTERISTICS  OF  NOZZLE 

In  nozzle  flow,  choose  the  incompressible  potential  flow  line 
and  the  corresponding  equipotential  line  as  coordinates  £,  n. 
Considering  that  the  flow  is  symmetrical  w.r.t.  the  central  flow 
line,  we  will  only  discuss  the  flow  in  the  upper  half  of  the 
nozzle.  Under  these  conditions,  the  boundary  conditions  for  the 
nozzle  can  be  expressed  as  follows: 

ij-oft  «Kf»o)  -  0  | 

i}  —  *?*)  — 

In  the  above,  (Ji^  is  the  amount  of  flow  passing  through  the  upper 

half  of  the  nozzle.  Because  of  the  limit  set  by  the  cross-section 

of  throat  in 4nozzle  flow,  cannot  be  any  arbitrary  value,  but  has 

a  maximum  value  which  we  call  , 

s 


% 
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For  the  sake  of  simplicity,  in  what  follows,  we  discuss  the 
case  where  the  cross-section  of  the  throat  is  also  symmetrical. 

We  will  see  later  that  the  results  can  be  extended  to  the  case 
of  unsymmetri cal  throat  cross-section.  In  the  nozzle  flew  cor¬ 
responding  to  symmetrical  throat  cross-section,  the  coordinate  £ 
for  the  cross-section  of  the  throat  is  0.  From  boundary  condi¬ 
tions,  Equation  (1G)  and  the  solution  for  the  stream  function 
tp(£,n),  Equation  (13),  one  obtains  for  the  amount  of  flow  passing 
through  the  cross-section  of  the  throat,  ^ ,  the  expression 


r  +  l  HKo,n) J 


(17) 


In  Equation  (17),  after  the  incompressible  flow  is  specified, 
H(0,n)  is  a  definite  known  function.  Therefore,  becomes  a  func¬ 
tion  of  g(0)  only.  When  ipfe  is  at  its  maximum,  one  has  [<ty>»/^K0)  1  “  0. 
We,  therefore,  differentiate  Equation  (17)  with  respect  to  g(.0)  to 
get 


]  [i 

J  l  r  +  1  v 


08) 


In  Equation  (18),  g(0)  is  multi-valued;  g  (0)  represents  the 
smallest  positive  real  root  and  is  the  flow-related  function  cor¬ 
responding  to  the  maximum  flow  ip  passing  through  the  cross- 

s 

section  of  the  throat.  Substituting  g  (0)  into  Equation  (17),  one 

s 

obtains  an  expression  for  the  maximum  flow  ip  as 

s 


T-l  rf(0)  1 

r  +  l 


i 


(19) 


In  the  nozzle  flow,  when  tp^  Is  smaller  than  ipg,  one  is  deal¬ 
ing  with  subsonic  nozzle  flow,  i.e.,  except  for  the  local  super¬ 
sonic  flow  that  might  appear  near  the  throat  wall,  the  flow  on  both 
sides  of  the  throat  is  subsonic.  When  ip^  equals  ip  t  one  has  super¬ 
sonic  nozzle  flow.  In  other  words,  under  this  condition,  the  flow 


15 
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in  the  throat  region  turns  from  subsonic  into  supersonic.  The 
attainment  of  maximum  amount  of  flow  in  the  nozzle  is  a  neccss- 

o 

ary  condition  for  guaranteeing  continuity  of  flow  in  the  throat 
region  and  avoiding  singular  points.  The  determination  of  max¬ 
imum  flow  tjj  is  thus  a  very  important  problem. 

b 


Let  the  boundary  conditions  be  such  that  ^  •  After 

is  given,  one  can  use  the  boundary  conditions  to  obtain  from 
Equation  (13)  the  integral  equation  for  determining  the  flow- 
related  function  g(c) 


■r*®  £*[•- 


XJZ±  _£££!_] 7-1 
r  +  1  //‘(f.OJ 


(20) 


This  is  an  integral  equation  for  g(.£)  with  £;  as  a  parametric 
variable.  For  a  given  E, ,  g(£)  in  Equation  (.20)  is  multi-valued. 

To  ensure  continuity  of  solution  and  that  the  density  remains  posi¬ 
tive,  the  value  for  g(£)  should  be  chosen  in  the  following  manner: 
For  subsonic  flow,  ^  <  tps,  take  g(£)  to  be  the  smallest  positive 
real  root  for  the  entire  range  of  £,  For  supersonic  nozzle  flow, 

^b  =  ^s*  the  choice  s(?)  is  determined  by  the  value  of  £ .  When 
C<0,  take  g(£)  to  be  the  smallest  positive  real  root;  this  corres¬ 
ponds  to  the  flow  in  the  contraction  section.  When  £  =  0,  take 
g(C)  to  be  the  smallest  positive  real  multiple  root;  this  corres¬ 
ponds  to  the  flow  in  the  cross-section  of  the  throat.  £>0,  take 
g(S)  to  be  the  second  real  root;  this  corresponds  to  the  flow  in  the 
divergent  section. 

When  y  =  1  +  and  v  is  a  positive  whole  number,  as  ^(£,q) 

can  be  expressed  as  a  polynomial,  the  corresponding  maximum  flow 

tj>  and  the  flow-related  function  g(£)  can  be  calculated  using  poly- 
s 

nomials.  Thus,  the  amount  of  flow  passing  through  the  cross-section 
of  the  throat  ( f;  =  0)  can  be  expressed  as 

* 

O.IJ.V^'CO)  (21) 

*•0 
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Differentiate  Equation  (21)  w.r.t.  g(0)  and  set  the  result  to 
zero.  One  obtains  the  algebraic  equation  satisfied  by  g(0)  at 
maximum  flow: 

2(2*+  l)/f*E*(0,.j»V*(0)  -0  (22) 

*-0 

Let  g  (0)  represent  the  smallest  positive  real  root  at  g(0)  in 
s 

Equation  (22)  and  substitute  it  into  Equation  (21)  to  obtain  bs 

9 

*.  “  S  VtW+'CO)  .  (23) 

k»  0 


The  relation  used  to  determine  g(.£)  from  boundary  conditions 
is  given  by 

9 

2  ,  t?*V*+,C£)  -  to*  (24) 

*-0 

In  Equation  (24),  g(C)  is  multi-valued.  The  choice  of  g(.£) 
under  various  conditions  is  the  same  as  in  Equation  (20). 

The  above  polynomial  form  for  the  stream  function  ip  ( .^ , n )  is 
applicable  for  the  cases  where  y  is  2,  1.5,  1.33,  1-25,  etc. 
Although  the  choice  of  y  is  thus  limited,  the  method  of  solution 
is  simple  and  easy  to  apply.  For  example,  when  y  is  2,  one  obtains 
from  Equations  (14),  (16),  (22),  (23),  (24),  the  following  analy¬ 
tical  expressions  for  g(£),  g  (0)  and  to  : 

s  s 

Kf)“  [2£/(0)/m(£)]  cos-j-  {st± cos-,[*u>(f)J }, 

«<?)  -  [E,(?,»fc)/E.(0,ij»)J* 

f/°)“  [Vt/E,(0,  »;*)]*  ’ 

<t>>  —  >?»*«(o),  *  -  to»/to. 


In  conformal  curvilinear  coordinates,  because  the  Lam£  coeffi¬ 
cients  are  equal,  the  relation  for  the  characteristic  line  is  the 
same  as  that  in  orthogonal  linear  coordinates.  Thus  we  have 

(*±#0  — tg "'(to«/to,),  M  “  *Sn“'  —  (26) 
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Usually,  the  characteristic  lines  obtained  by  taking  the 
positive  sign  in  Equation  (26)  are  called  the  first  family  of 
characteristic  lines,  and  those  obtained  by  taking  the  negative 
sign  are  said  to  belong  to  the  second  family. 

For  convenience,  the  second  family  characteristic  line  that 
passes  through  the  point  (£  —  —  on  the  cross-section  of  the 

throat  is  called  the  influence  line.  From  the  characteristics  of 
supersonic  flow,  the  flow  upstream  of  the  influence  line  will  not 
be  affected  by  the  wall  surface  of  the  divergent  section. 

The  second  family  characteristic  line  that  passes  through  the 
sound  velocity  point  on  the  center  line  (f  —  tj  —  o)  is  called  the 
limiting  characteristic  line.  The  flow  downstream  the  limiting  w 
characteristic  line  is  generally  considered  not  to  affect  the 
shape  of  the  sound  velocity  line.  In  nozzle  flows,  the  wall  point 
of  the  limiting  characteristic  line  is  upstream  the  cross-section 
of  the  throat,  and  the  velocity  line  is  affected  by  the  maximum  flow 
,  which  is  dependent  on  the  shape  of  the  throat  wall  at  the  cross- 

O 

section.  Hence,  when  adjusting  the  wall  shape  of  the  throat  down¬ 
stream  the  limiting  characteristic  line,  one  should  keep  unchanged 
so  as  to  ensure  that  the  shape  of  the  velocity  line  stays  the  same. 

Usually,  the  first  family  characteristic  line  that  passes 
through  the  sound  velocity  point  on  the  center  line  is  called  the 
branchline.  When  using  the  hodograph  method  to  discuss  nozzle 
flow,  the  branchline  is  a  singular  line  of  the  hodograph.  It  sat¬ 
isfies  the  relation  I0(?,  0)/d(tf>,  </>)]  —  0lr). 

♦ 

To  find  the  time  x(s,i^)  spent  by  a  gas  mass  point  on  any 
stream  line,  let  s  =  Sq  be  the  zero  reference  time,  i.e.,  tCsq,i|i) 

=  0.  Then  the  time  spent  along  any  stream  line  ip  =  C  can  be 
written  as 

rO,  -  ['  -~r,  f(S,  O  -  f*  HO  +  4>\/ 4>X)* d£  <27> 
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in  which  £  is  the  value  of  £  on  each  streamline  corresponding  to 
Sq.  t(s  ,</>)  is  dependent  on  the  geometric  shape  of  the  nozzle.  In 
the  design  of  fast  acceleration  short  nozzles,  the  v/all  contours 
of  the  nozzle  selected  has  to  satisfy  the  requirement  of  minimum 
x(s  ,<//). 


V.  NOZZLE  COORDINATES 


According  to  the  theory  of  complex  variables  and  that  of  in¬ 
compressible  flow,  choose  symmetrical  hyperbolic  wing  shape  in  the 
incompressible  flow  velocity  plane  and  obtain  a  symmetrical  nozzle 
in  the  physical  z  plane.  The  transformation  relationship  equation 
f(0  between  z  plane  and  complex  potential  ;  plane  is  given  by 


yj 

chy(A  +  ol 

Z  "  —  < 

?  —  (1  —  <r)coch  A  In 

L 

- — — 

ch-i-U-f)  j 

(28) 


in  which  A, a,  are  shape  factors  of  the  nozzle,  and  n  is  the  con¬ 
traction  ratio  of  the  nozzle.  Details  will  be  given  in  another 
paper.  Separate  the  real  and  imaginary  parts  of  Equation  (28), 
and  obtain  the  following  relation  between  the  x,y  coordinates  and 
the  £,n  coordinates: 


*(S>  v) 

v(5,  >?) 


±{i 

; - — (1  —  g)cothAln 

[ch  (A  +  f)  4-  cos 

•J*  1 

2 

Ich  (A  —  f)  +  COS: 

~[v 

—  (1  —  cOcothAtg-1  ( 

sh  A  sin  t;  \1 

Vt  «■ 

*ch  £  +  ch  A  cos  tj/  J 

(29) 


From  Equations  (6d),  (28)  or  (29),  one  obtains  the  amplifi¬ 
cation  factor  H  ( £ ,  p )  of  z  plane  w.r.t.c  plane 


W(£» »?) 


ZL  [ (  cost;  4-  gch  Ach  H)3  —  (oJ  cha  A  —  1)  sh* 
t;*  L  (cost;  +  chi ch^)’  —  sh1 Ash' {  J 


(30) 


When  p  =  p^,  Equation  (29)  gives  the  contour  line  of  the  wall 
of  a  symmetrical  nozzle,  with  £  as  a  parametric  variable.  From 
analytical  geometry,  the  following  relation  exists  among  the 
nozzle  wall  shape  parameters  r  (radius  of  curvature  of  throat  wall). 
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n  (contraction  ratio),  0^  (largest  an  g  1  e  of  inclination  of  wal] ) 

and  f,  (point  of  inflection  at  the  wall): 

± 


r  —  [(1  +  y')W/y")^o.  n  —  y(oo,  TJ,) 

0/  “  'g-yCfi.  »?*),  y"(ft»  *2*)  —  o 


OD 


In  the  above  equations,  y*  and  y"  represent  the  first  and  second 
derivative  respectively  of  y  w.r.t.  x.  From  Equations  (30),  (3D, 
and  n  =  ,  we  obtain  the  following  relations  between  the  wall 

shape  parameters  and  X,o,  : 


n  (cost ]b  -F  gch  1)* 
ij»  (l  — c)chlsin»;» 


n  “  [  1  —  - - -  COth  1  tg 

1  •?» 

■[ 


e,  —  tg 


Utff,  +  0  +  <0 


r-i  (  shlsin  i?» _ \1~‘ 

M  *4*  chi  cosij,/J 

(1  —  <r)  ch  1  sin  sh  g, _ 1 

■  g)  ch  icosij,  ch  ft  +  <rchJ  1  —  sinJ»j,J 


ch  ft  -  [A  C«h»l  +  cos’,*  +  1)1  f  X  cos-l{cos-»  fCl  ,t  *>.«=, h  1 

14  3  [A(achU  +  co^  +  l)|fJ 


(32) 


With  given  nozzle  shape  parameters  r,  n,  6^,  one  can  obtain 
from  the  above  equation  the  corresponding  value  for  A,o,n^  and 
hence,  the  desired  nozzle  shape  line. 

VI.  EkU,n),  FUNCTION  RELATING  NOZZLE  TRANSFORMATION 

For  the  above  nozzle  which  had  adjustable  radius  of  curvature 
r,  contraction  ratio  n  and  largest  angle  of  inclination  of  wall  0^ , 
when  y  =  1  +  where  v  is  a  positive  whole  number,  the  corres¬ 
ponding  function  relating  to  the  transformation,  Ek(£,n)  can  be 
obtained  in  the  form  of  an  analytical  expression  by  integration. 
Substitute  Equation  (30)  into  Equation  (Db),  express  the  inte¬ 
grand  in  terms  of  partial  fractions,  integrate  separately  and  sum 
together  to  obtain  Ek(g,  n): 


ft 
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E«d,  v) 


Cod,  v)  *“  *1 

G,(£,  *;)—  2  [Si«d»  a)‘V»d,  *?>  «)  +  B'»d.  -“WmiS,  r],  -a)]  , 

»■  1 

cha  *“  achl 


(33) 


In  Equation  (33),  Nm(£,n,ot)  and  Blm(£,a)  are  respectively: 


W.d,  «?,  a)  -  77"  'K-1  (th  *  r6  -7) 

thf  4*  tha  V  2  2  / 


N«(|,»/, «)  “ 


(m  4-  1  )(th  |  4-  th  a 
(m  —  17,  a)  — 


y  (2m  -  3)(1  +  th  a)  th  17,  a) 


tech  a  tech  j;  sin  t; 


(I  4-  th  a  th  £  4-  tech  a  tcch  | cos  17)"*  1 
Bud,  «)  *■  1  4*  coth  a coth  £  ^rhJ£  —  -  -  -  -  tech1  Sj 

*-«. «)  -  E  c-ir-  ^  ocssr* 

/,  m  ■»  1 ,  2,  3,  •••;«<!/ 


J 


(34) 


For  convenience,  the  expressions  for  B-Lm(C,ot)  for  1  ,m  $  3 
are  found  from  Equation  (3^d)  and  listed  below.  These  correspond 
to  the  cases  where  y  =  1.5  and  1.33: 


S ad*  o)  —  —coth  a  coth  $B„d,  a)Bu(£,  —  a) 

B„(f,  <0  -  S?,d,  «) 

®n(f,  <0  “  j  coth  a  coth  £B„(£ ,  a)[  B„d,  a)  —  fl„(f,  —  o)] 

S«d,  a)  -  i-  a)  B„d,  «) 

S„d,  a)  -  Bud,  <0 


(35) 
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At  the  cror,r.-3cct ion  of  the  throat,  C  =  0.  Use  this  r<  lot  ion 
in  Equations  (33)  and  3Jlb);  integrate  to  obtain  the  less  compli¬ 
cated  expression  for  E^(0,n): 


Ek(0,  ./)  -  (z*y  2  (I7^),'^(0,  ,),  *  -  0,  1,  2,  ... 

A’«(0,  v')  “  »?>  AA(0,  l)  “  2corh  a  tg-1  ^th  ■—  •  tg 


N-(0,  V) - 


ch"-!a 


coth1  a  j  (2m  —  B)A’m_,(0,  »;) 

seeh  a 


m  —  1 

-Cm-  2)A/m_,(0,  ,)- 
*n  “*  2,  3,  •  •  •  \  m  ^  ^ 


(1  +  sech 


■linn  _] 

acosij)”-'  J 


Li) 


Usually,  ochA>l,  and  cha  is  greater  than  1,  where  a  is  a  real 
number.  When  the  radius  of  curvature  of  the  throat  wall  is  very 
small,  for  instance,  when  nb  =  ir/2,  r  =  0.3,  achX  <1.  In  this 
case,  a  is  an  imaginary  number.  The  corresponding  N  (.£,n,u)  and 
Bim(C,a)  are  complex  numbers.  Considering  that  H(.^,ri)  is  always  a 
positive  real  number,  E^(£;,n)  must  also  be  a  positive  real  number. 
When  a  is  imaginary,  let  a  =  iB,  in  which  B  is  a  real  number.  Then 
one  has  cha  =  cos  B,  sha  =  i  sin  B,  and  the  identity 


+ 


1  ±  ,’f  .  re  1\  =  JL  tK-‘  f _ ±lsAA _ 

2  2  /  2  'cos»j  +  cos/J  tech  f 

JL  th-i  (  —JiEJLlBl _ ) 

2  \  cos  tj  +  sec  £  ch  §  / 


) 


(37) 


Thus,  using  the  above  identify  and  equations  to  separate 
N-^(C,n>a)  and  B^(£,a)  into  real  and  imaginary  parts  and  writing 
Nm(C,n,a)  and  Blm(C,a)  as  functions  of  11^(5, n, a)  and 
after  the  results  for  Nm(5,n,a)  and  (C,a)  are  substituted  into 
Equation  (33),  one  obtains  the  corresponding  G^(£,n)  and  E^(C,n). 
Furthermore,  when  ochA  <1,  one  can  use  the  integral  Equations  (13) 
and  (17)-(20)  to  obtain  the  results  directly. 
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VII.  AH 


OF  MCZZLK  FLOW  FIELD  CO*':  STATION 


In  this  section,  wc  show  computations  based  on  the  general 
solution  for  rissole  flow  given  above  for  the  simplest  case,  i.e., 
when  y  equals  2.  From  Equations  (15)  and  (25),  one  obtains  for 
Y  =  2  the  simple  analytical  expression  for  tl(S,n): 


I/)  -  [  3  w(£)  -  E,Q,  v)gX 4)1/2  > 

e(S)  -  {!r±coS-,[^£r1/!C0,  ,?»)£««(*,  ^)]} 


v  +  (^)  y,  a) 

+  Bn(£,  —  n>  —a)] 

Bn(? ,  a)  “*  1  +  coth  a  coth  |^thJ  £  —  -  ~  scch1 

~  u  tg'1  (rh  £±«  .  rg  A) 

th{  +  th«  \  2  2/ 

0.  *“  y*),  %  “  < P>/4>„  ch  a  ■»  a  ch  X 


(3b) 


Under  these  conditions,  p(£;,n)  in  Equation  (7a)  can  be  ex¬ 
pressed  as 

pc r,  o  -  r ...  {a  co.-  f  i  -  u\  +  *»]}  (39) 

In  Equation  (39),  c°*  1  [ 1  ^  (</>{  4-  <^5)j  has  yalues  between  0  and 

it  in  the  subsonic  range,  and  between  it  and  2tt  in  the  supersonic 
range . 

The  flow  for  the  case  in  which  y  =  2  is  equivalent  to  the 
motion  of  gas  which,  in  the  water-table  experiment,  was  simulated 
with  the  free  motion  on  the  water  surface  [933-  For  the  cases 
where  y  equals  1.5,  1.33,  1.25,  1.20,...,  the  stream  function  f'(t;,n) 
has  analytical  expressions  similar  to  Equation  (38);  for  any  other 
values  of  y,  one  can  use  the  integral  expressions  Equations  (13) 
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the  r^'-sul tr>  for  Copl¬ 


and  (20).  As  wo  aro  1  ini  tod  in  space  hoz>>-, 
putation  for  other  values  of  y  will  bo  given  olscwboro. 


Whc-n  cha  =  och.\<l,  the  corresponding  and  N  (^,n,oi) 

in  Equation  (38d,e)  are  complex  numbers.  In  this  case,  let  a  = 
i3,  use  Equation  (37)  to  separate  B^(f.  ,a)  and  N^(C,n,a)  into  real 
and  imaginary  parts,  and  substitute  into  Equation  (38c).  One  then 
obtains  the  corresponding  expression  for  ( f.  ,n): 


E.a,*?)  -  fey  (tg1  fi  th‘  O  -  thJ  i) 

X  coth  f  rg-»  ( - *12*1 - )  +  (2tcc>fi-l±£  Kch>  s  )  agj} 

Vcosjj  4-  cos/J  sech£/  V  a  k 


th-(. 


X  th* 


sin  rtgft 


cos  i;  +  sec  fich  ^ 


;)M?)' 


(40) 


One  can  use  Equations  (38),  (39),  (7b),  (4b)  and  Equations 
(26)  and  (27)  to  obtain  the  density  p(£>n)  in  the  flow  field, 
velocities  u(£,n),  v(£,n)  and  Mach  number  M(.£,n)  distribution, 
along  with  the  influence  line,  limiting  characteristics  line, 
branchline  and  equitemporal  lines. 


Note  that  the  stream  function  given  in  Equation  (38)  contains 
three  shape  factors  l,o,nb  and  one  flow  factor  k.  For  given 
nozzle  wall  shape  yarameters,  r,  n,  and  6-^,  one  can  determine  the 
corresponding  shape  factors  A,a,n^  from  Equation  (.32-)  <  Flow  fac¬ 
tor  k  varies  in  the  range  between  0  and  1.  k  less  than  1  gives  sub¬ 
sonic  nozzle  flow,  while  k  equal  to  1  gives  supersonic  nozzle  flow. 


Computation  for  a  typical  nozzle  is  given  below.  The  shape 
parameters  of  the  nozzle  are  taken  to  be,  respectively,  r=2,  n=2, 
61=23.3°;  the  corresponding  shape  factors  are,  respectively,  \=2,23, 
<7=0.434  ,  nb=  7j-.  The  calculations  are  done  for  k  equal  to  0.50,  0.90, 
0.98  and  1.00,  respectively.  These  correspond  to  lower  subsonic, 
medium  subsonic,  higher  subsonic  and  supersonic  nozzle  flow,  res¬ 
pectively.  The  results  of  the  calculations  are  shown  in  Figures 
1-6. 
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Figure  1  riv"'3  the  distribution  of  linos  of  equal  Mach  number 
for  subsonic  nozzle  flow  (k=0.50).  The  largest  Mach  number  in  the 
figure  is  0.335.  Hence,  we  can  consider  the  flow  to  be  incompress¬ 
ible. 

Figure  2  gives  the  distribution  of  lines  of  equal  Mach  number 
for  medium  subsonic  nozzle  flow  (k=0.90).  It  can  be  seen  from  the 
figure  that  the  Mach  number  M>0 . 7  near  the  throat  wall.  Hence,  one 
must  take  into  consideration  the  effect  of  compressibility  in  this 
region . 

Figure  3  gives  the  distribution  of  lines  of  equal  Mach  number 
for  higher  subsonic  nozzle  flow  (k=0.98).  It  can  be  seen  from  the 
figure  that  local  supersonic  flow  appears  near  the  throat  wall. 

This  type  of  flow  also  belongs  to  the  realm  of  transonic  flow 
problems . 

Figure  M  gives  the  distribution  of  lines  of  equal  Mach  number 
for  supersonic  nozzle  flow  (k=l).  It  also  shows  the  influence  line, 
limiting  characteristic  line  and  branchline.  AE  represents  the 
sound  velocity  line.  Note  that  the  sound  velocity  point  A  at  the 
wall  is  located  upstream  of  the  cross-section  of  the  throat  and  the 
sound  velocity  point  E  on  the  center  line  is  located  downstream  of 
the  cross-section  of  the  throat.  BE  is  the  limiting  characteristic 
line.  It  is  the  second  family  characteristic  line  that  passes 
through  point  E.  Wall  point  B  is  located  upstream  of  point  C  on 
the  cross-section  of  the  throat  and  downstream  sound  velocity  point 
A.  It  is  usually  regarded  that  the  flow  field  downstream  of  the 
limiting  characteristic  line  and  the  shape  of  the  wall  will  not 
affect  the  shape  of  the  sound  velocity  line.  * 

In  nozzle  flows,  the  entire  flow  is  affected  by  the  maximum 
flow  <Ji  ,  which  in  turn,  is  dependent  on  the  shape  of  the  throat 
wall.  Hence,  to  ensure  that  the  shape  of  the  sound  velocity  line 
will  not  vary,  one  must  make  sure  that  changing  the  wall  shape 
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»»2.0,  f— 2.0.  0,  — 23.3*.  r«2.0,  *-0.9« 


Figure  *1 .  Characteristic  lines  for  supersonic  noccio  flow 

—  2-0,  r=> 2.0,  S,=23.3*,  r=2.0,  4-J.O 


- lines  of  equal  Mach  number  sound  velocity  line 

- •-limiting  characteristic  line 


AE 


Figure  5.  Distribution  of  Mach  number  at  cross-section  of  throat 

—  2.0.  —2.0,  e,  =  23.3\  r-2.0 

1 — Mach  number;  2 — distance  from  center  line,  y 
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Figure  6.  Equi temporal  lines  of  supersonic  nozzle  flow 

■  -2.0,  r= 2.0.  0,  =  23.J\  r  =  2.0,  *=»t.O 

downstream  the  limiting  characteristic  line  v;ill  not  change  the 
maximum  flow  .  CF  is  the  influence  line,  i.e.,  the  second 

O 

family  characteristic  line  that  passes  through  wall  point  C  at 
the  throat.  The  flov;  downstream  CF  will  not  affect  the  flow  in 
the  convergent  section.  In  other  words ,  the  flow  upstream  CF  is 
not  affected  by  the  wall  of  the  divergent  section.  Therefore,  one 
can  choose  a  line  downstream  of  the  sound  velocity  line  and  up¬ 
stream  of  the  influence  line  as  the  initial  condition  for  the 
design  of  the  divergent  section.  According  to  the  requirements 
of  the  initial  value  problem,  this  line  cannot  be  a  characteristic 
line.  ED  is  the  branchline,  the  first  family  characteristic  line 
that  passes  through  the  sound  velocity  point  on  the  center  line. 

When  using  the  hodograph  method  to  analyze  nozzle  flow,  the  branch¬ 
line  is  a  singular  line  on  the  hodograph  plane  [27]. 

Figure  5  gives  the  distribution  of  Mach  number  in  the  cross- 
section  of  the  throat  for  different  amounts  of  flow.  From  the 
figure,  one  can  see  that  the  greater  the  flow  factor  k,  the  less 
uniform  the  distribution  of  Mach  number.  Hence,  the  assumption 
of  a  linear  sound  velocity  line  in  design  of  supersonic  nozzles  £83] 
is  not  reasonable.  One  must  take  into  consideration  the  effect  of 
a  curvilinear  sound  velocity  line,  especially  for  the  case  of 
smaller  radiu's  of  curvature  of  throat  wall. 
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Figure  C  gives  tho  equiti  inporal  lines  of  supersonic  nozzle 
flow.  The  zero  reference  time  is  taken  at  £  =  f,^.  From  the 
figure,  one  can  see  that  the  density  of  equi temporal  lin-s  is 
higher  in  the  subsonic  region,  i.e.,  the  gas  spends  more  time  in 
the  subsonic  region.  Hence,  to  hasten  the  flow  of  gas  in  the 
nozzle  so  as  to  increase  the  freezing  efficiency  of  high-temperature 
gases  under  non-equilibrium  conditions,  one  must  give  heed  in  the 
design  of  the  contour  lines  of  the  wall  in  the  convergent  section. 

VIII.  Discussions 

1)  The  solution  for  subsonic,  transonic  and  supersonic  nozzle 
flow  has  been  given  in  a  generalized  form.  It  is  applicable  to 
lower  subsonic,  medium  subsonic,  higher  subsonic  and  supersonic 
nozzle  flows. 

2)  We  have  given  a  nozzle  whose  wall  shape  can  be  chosen  as 
desired  and  its  solution.  The  wall  shape  of  the  nozzle  Is  deter¬ 
mined  by  three  parameters:  radius  of  curvature  r  of  the  throat 
wall,  contraction  ratio  n  and  largest  angle  of  Inclination  01  of 
wall,  r  can  be  taken  to  be  any  value  from  0.05  to  infinity  and  n 
can  have  any  value  from  1  to  infinity.  In  choosing  0^,  the  upper 
limit  increases  with  increasing  contraction  ratio  n,  and  increases 
with  decreasing  radius  of  curvature  r  of  throat  wall.  For  example, 
when  r=l,  and  n  is  2,5,10,20,  the  upper  limit  for  0^  is  res¬ 
pectively  23°,  ^5°,  57°,  60.5°. 

3)  The  solution  for  nozzle  flow  given  is  valid  for  any  value 

of  the  ratio  of  specific  beats  y.  Under  general  conditions,  the 
solution  for  the  stream  function  is  in  integral  form.  When 

y=  1  +  (1/v),  where  v  is  a  positive  whole  number,  ip ( C » n )  takes  the 
form  of  an  odd-powered  polynomial. 

i| )  The  expression  for  the  maximum  flow  ^  has  been  given,  ^ 

s  s 

is  a  function  of  the  shape  of  the  wall.  When  the  radius  of  curvature 
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r  of  4  he  1  hrr<-:t.  wall  is  largo,  v  ^  1.  With  decreasing  r, 

s  ’  s 

decreases . 

5)  The  results  of  calculation  show  that  -when  r  is  small,  the 
sound  velocity  line  is  far  from  being  linear.  Under  these  condi¬ 
tions,  it  is  not  reasonable  to  use  a  linear  line  for  the  sound 
velocity  line  to  be  used  as  initial  condition  in  the  design  of 
the  divergent  section  [82,83]-  A  similar  situation  prevails  in 
axially  symmetric  nozzle  flows  [75,76]. 

6)  The  second  family  characteristic  line  that  passes  through 
the  wall  point  (£  =  0,  n=  n^)  in  the  cross-section  of  the  wall  is 
called  the  influence  line.  From  the  properties  of  supersonic  flows, 
one  knows  that  the  flow  upstream  of  the  influence  line  is  not 
affected  by  the  shape  of  the  wall  in  the  divergent  section.  Hence, 
one  can  take  a  line  downstream  of  the  sound  velocity  line  and  up¬ 
stream  of  the  influence  line,  specify  the  magnitude  and  direction 

of  the  velocity  on  this  line,  and  use  it  as  initial  condition  in 
the  design  of  the  divergent  section.  By  the  requirements  of  the 
initial  value  problem,  this  line  cannot  be  one  of  the  characteris¬ 
tic  lines. 

7)  The  solution  given  for  the  simplified  equation  for  the 
stream  function  has  a  high  degree  of  accuracy  in  the  transonic 
region.  It  has  also  given  satisfactory  results  for  the  higher  sub¬ 
sonic  region  and  lower  supersonic  region.  For  higher-order  solu¬ 
tions,  one  can  take  this  solution  as  a  basis  and  use  the  general 
form.  Equation  (.8),  for  the  stream  function  in  conformal  curvi¬ 
linear  coordinates  to  dq  the  analysis, 

8)  The  solution  given  for  the  subsonic  nozzle  flow  corres¬ 
ponds  to  that  for  symmetrical  throat  cross-section.  One  can  further 
discuss  higher  order  approximations,  using  this  solution  as  a  basis. 

9)  For  nozzles  with  hyperbolic  or  circulnr-arc-shaped  throat 
wall,  the  solution  given  for  the  transonic  flow  has  a  simple  form 
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and  is  convenient  to  use  in  t h"  analysis  of  the  pronert  cif 
flow  in  the  throat  region. 

10)  The  characteristics  of  transonic  flow  in  a  hyperbola¬ 
shaped  nozzle  obtained  by  our  method  are  in  good  agreement  with 
those  obtained  by  Cherry  [52]  and  Serra  [8l]  for  the  case  where 
y  =  1.4.  The  results  will  be  given  elsewhere. 

11)  From  Equation  (20),  it  is  obvious  that  when  q  is  close 

to  1,  y  has  little  effect  on  g(£).  Roughly  speaking,  in  this  case, 
g(5)  is  dependent  on  the  geometric  shape  of  the  nozzle  only, 

Hence,  for  a  given  nozsle,  the  shape  of  the  sound  velocity  line  is 
basically  unaffected  by  the  ratio  of  specific  heats  y. 

12)  To  study  the  effect  of  nozzle  wall  shape  on  the  flow  field, 
particularly  the  relation  between  the  flow  characteristics  of  the 
transonic  region  at  the  throat  and  the  geometric  shape  of  the 
nozzle,  one  can  use  the  nozzle  given  in  this  paper,  do  calculations 
using  various  shape  parameters  r^S-^,  and  make  a  table  for  refer¬ 
ence  in  the  engineering  design  of  nozzles. 

13)  When  it  is  not  feasible  to  express  in  simple  analytical 
forms  the  incompressible  potential  flow  of  the  nozzle,  one  can  use 
numerical  methods  to  calculate  the  streamlines  and  equipotential 
lines  of  the  incompressible  flow.  On  this  basis,  one  can  use  the 
present  method  to  compute  the  corresponding  subsonic,  transonic 
and  supersonic  nozzle  flow. 

14)  This  method  can  be  extended  to  flow  around  objects.  Here, 
the  manner  in  which  the  flow-related  function  is  determined  from 
the  boundary  conditions  differs  from  that  in  the  case  of  internal 
flows.  We  will  discuss  this  in  another  paper. 

15)  In  analyzing  axially  symmetric  flows,  one  can  use  the 
streamlines  and  equipotential  lines  of  the  corresponding  incom¬ 
pressible  potential  flow  as  coordinates  to  solve  for  the  axially 
symmetric  flow  of  compressible  gases. 


tx.  oorjoujnio:.':'. 

Tlii".  art i  e ]  *.•  r;1vf'$  th->  general  Corn  of  the  equation  for  the 
two-dimensional  constant  steady  irrotattonal  isentropic  flow  of 
an  ideal  gas  in  conformal  curvilinear  coordinates.  Simplified 
equation  for  the  stream  function  and  its  general  solution  is  given 
using  the  streamlines  and  cquipotent ial  lines  of  the  corresponding 
incompressible  potential  flow  as  coordinates. 

The  results  are  applied  to  nozzle  flow,  giving  the  solution 
for  nozzle  flow  from  subsonic  through  transonic  to  the  supersonic 
region,  for  any  desired  choice  of  the  wall  shape,  i.e.,  radius  of 
curvature  of  throat  wall,  contraction  ratio  and  largest  angle  of 
inclination  of  wall.  This  solution  is  valid  for  different  ratios 
of  specific  heats.  As  an  example,  we  have  calculated  the  flow  char¬ 
acteristics  of  a  typical  nozzle. 

In  this  paper,  the  relation  between  the  maximum  flow  and  the 
shape  of  the  wall  is  given.  The  maximum  flow  decreases  with 
decreasing  radius  of  curvature  of  throat  wall.  For  very  large 
radius  of  curvature  of  throat  wall,  e.g.,  r>5,  the  maximum  flow 
can  be  approximated  by  1. 

This  paper  gives  the  flow  characteristics  of  a  nozzle  for 
different  amounts  of  flow.  The  case  with  flow  factor  k  less  than 
1  corresponds  to  subsonic  nozzle  flow,  i.e.,  except  for  local  super¬ 
sonic  region  that  may  appear  near  the  throat  wall,  the  flow  on  both 
sides  of  the  throat  is  subsonic.  When  k  equals  1,  the  flow  is 
supersonic,  i.e.,  near  the  throat,  the  flow  changes  from  subsonic 
to  supersonic.  * 

In  supersonic  nozzle  flows,  the  shape  of  the  sound  velocity 
line  is  greatly  affected  by  the  contour  of  the  wall,  especially 
the  radius  of  curvature  of  throat  wall.  For  a  nozzle  of  given  wall 
contours,  the  shape  of  the  sound  velocity  line  is  not  much  affected 
by  the  ratio  of  specific  heats,  y.  The  results  of  calculations 
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for  hyperbola-shaped  nor/.:]  os  show  that  v.'hen  y  varies  fror:  1.2 
to  1.667,  the  sound  velocity  point  on  the  center  line  moves  down¬ 
stream  a  little  bit,  the  sound  velocity  point  on  the  wall  moves 
upstream  a  little  bit,  but  the  changes  are  very  small.  For  details 
see  another  paper  to  be  given  later. 

The  second  family  characteristic  line  that  passes  through  the 
wall  point  of  the  cross-section  of  throat  is  called  the  influence 
line.  Considering  that  the  flow  upstream  of  the  influence  line  is 
not  very  much  affected  by  the  wall  contours  of  the  divergent  sec¬ 
tions,  one  can  choose  a  line  between  the  sound  velocity  line  and 
the  influence  line,  use  the  method  given  in  this  paper  to  give  the 
magnitude  and  direction  of  velocity  on  this  line  and  use  this  as 
initial  condition  for  the  design  of  the  divergent  section.  Down¬ 
stream  the  influence  line,  the  contour  line  of  the  divergent  sec¬ 
tion  can  be  designed  by  means  of  the  method  of  characteristic  lines 
dy  the  requirements  of  the  initial  value  problem,  this  initial  line 
cannot  be  a  characteristic  line. 

This  method  can  be  extended  to  flow  around  an  obstacle  oyer 
weir  flow  in  the  nozzle  and  flow  through  a  grating.  Particularly, 
one  can  obtain  better  results  in  the  transonic  region.  Furthermore 
the  method  can  be  extended  to  cases  where  there  are  chemical  reac¬ 
tions  under  equilibrium  or  non-equilibrium  conditions,  or  where 
there  exists  axial  symmetry. 
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RESEARCH  ON  THE  HISTORY  OF  MECHANICS 


Normal  College  of  Inner  Mongolia  '  Li  Ti 

Through  production  struggle,  class  struggle  and  scientific 
experimentation,  the  working  people  of  our  country  in  ancient 
times  made  great  contributions  to  mankind  v/ith  their  achievements 
in  science.  Our  people  are  good  at  utilizing  accumulated  exper¬ 
ience  and  skillfully  making  various  tools  and  machines  based  on 
mechanical  principles.  They  have  invented  many  things,  one  of 
which  is  the  rocket. 

A 

I.  THE  ANCIENT  "ROCKET" 

"Rocket"  is  a  term  that  appeared  early  in  our  country.  In  the 
concise  edition  if  Wei  Liieh ,  it  is  recorded  that  when  Chu  Ke  Liang 
was  attacking  He  Chao,  he  "raised  the  tall  ladders  and  had  his 
men  climb  onto  the  wall.  Chao  used  rockets  to  shoot  at  the  ladder. 
The  ladder  caught  fire  and  all  the  men  on  it  were  burned  to  death". 
Tu  Huey  Tu  Chuan,  vol .  92  of  Sung  Shu,  records  the  battle  between 
Tu  Huey  Tu  and  Jui  Hsiin  on  water,  in  which  "Huey  Tu  personally 
fought  on  the  battleship,  had  rockets  shot  onto  Hsiin's  battleships 
which  all  caught  fire.  Hsiin  lost  the  battle  immediately.  He  was 
shot  by  an  arrow  and  died  in  the  water".  In  Wang  Szu  Chonr  Chuan, 
vol.  62  of  Pei  Shj_h ,  there  is  also  an  account  of  how  "Szu  Cheng 
used  rockets  to  shoot  "the  attacking  enemies  and  burned  their 
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attackin'*  i  tia'  :  ■>  i  te  a  "  .  To-  :’<■  a  er-  u;g  are  •  •videnro  that  th'  f'- 
existed  hln  '  call  'd  " cia-t  "  as  ..-arly  ns  In  lh<*  -  ra  of  t  ho 

Three  v.a:  -b  jont,  Inu*  <1  t  o  be  used  during  the-  Ooui-h  and 

North  l)y* T  ••/  .  7h'»:*--  "r~  rlo'S  s- "  v:*1 r>'*  rh^t  out .  Apparent  ly  , 

those  were  just  arr-  ws  with  f  1  .q  mm.  able  substances  addl'd  to  burn 
the  object  a.  an  ■  n  :  no’  the  kind  of'  rocket  that  is  propelled  by 
explosives . 

This  type  of  "rocket"  was  further  developed  in  the  Tang  Dynas¬ 
ty.  In  Vol .  ^  of  "Chen  Chi  Chih  Ti  Tai  Pai  Yin  Ching",  Li  Ch'uan 
recorded  a  kind  of  rocket  in  greater  detail.  "Rockets:  Attach  a 
ladle  of  oil  to  the  arrowhead  and  shoot  the  arrow  at  the  fortress. 
When  the  ladle  breaks  and  the  oil  leaks  out,  follow  up  with  a 
rocket  so  that  a  fire  is  started.  Then  follow  up  with  more  oil- 
'  ladles  and  the  fortress  is  all  burned  down".  This  account  may 
have  been  taken  from  the  Ch'ing  Dynasty  edition  of  Wei  Kung  Ping  Fa 
edited  by  Li  Ching  during  the  early  period  of  the  Tang  Dynasty. 

The  wording  is  almost  identical  to  that  used  to  describe  "rockets" 
in  Tai  Pai  Lin  Ching  This  kind  of  rocket  is  nothing  but  an 

ignited  arrow.  A  fire  starter  is  shot  at  the  object  to  ignite  it. 
To  make  a  big  fire,  combustible  oily  substance  precedes  the  rocket, 
and  after  the  fire  is  started,  more  oil-ladles  are  added  to  feed 
and  expand  the  fire. 

In  the  Tang  Dynasty  there  was  another  kind  of  fire  arrow  which 
employed  a  bow  to  shoot  out  the  arrows  that  carried  fire.  Details 
were  also  recorded  in  Tai  Pai  Yin  Chin,  vol.  4 .  It  goes  as  follows; 
"Fire  arrow :  Those  who  can  shoot  300  feet  away  attach  a  ladle  of 
fire  onto  their  arrowhead  and  several  hundreds  of  such  arrows  are 
shot  into  the  enemy's  tents.  Their  supplies  will  be  burned  away. 

One  then  takes  advantage?  of  the  confusion  to  launch  an  attack". 

Thus,  a  ladle  of  fire  is  attached  to  the  arrowhead  and  used  as  a 


Wo  1 _ Kung  _Pi  n<~  F-q  edited  by  Li  Ching  is  lost.  So  vol.  3  of  We  1 

Kung  Plug  i''a  Chi  ion,  edited  by  Wang,  Tsung  I  (Ch'ing  Dynasty)’." 
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fire  starter  to  burn  Mi'  object.  Thin  In  different  from  the 
above-mentioned  rocket,  which  first  thrusts  off  oil  indies  to  aid 
combustion.  However ,  these  are  similar  in  principle  in  that  a 
fire  starter  is  slur."  cut  via  elastic  force. 

Since  the  Sung  Dynasty,  rockets  have  found  more  and  more  fre¬ 
quent  use  on  battlefields.  "Ping  Kao"  vol.  131  of  Hsu  l>n  Hsien 
T'ung  Kao  recorded  that  in  April  of  the  8th  year  of  Tien  Huey's 
rule  of  Chin,  (1130  A.D.),  there  was  a  battle  between  Chin  and  Nan 
Sung.  The  Chin  soldiers  "shot  rockets  from  the  light  boats.  Smoke 
and  flame  covered  the  whole  place  and  the  enemy  was  badly  defeated". 
The  following  year,  in  the  battle  of  Tang  T'u  between  Chin  and 
Sung,  the  Chin  soldiers  "shot  rockets  and  burned  the  lookout  tower" 
of  Sung's  city  wall.  (Yiieh  Yen  Jo:  Yiin  Lu  Man  Ch'ao,  vol.  7). 
Similar  accounts  abound.  For  example,  in  Hsianm  Hsien  Shou  Cb'eng 
Lu,  it  is  recorded  that  in  April  of  the  2nd  year  of  Kai  Hsi  (1206 
A.D.),  the  army  defending  Hsiang  Hsien  used  "explosive  arrows  to 
shoot  and  burn"  the  Chin  soldiers.  And  in  "Ping  Kao"  vol.  131  of 
Hsu  Wen  Hsien  T'ung  Kao,  there  is  an  account  of  the  "fire  arrows" 
used  by  the  Mongolian  army  in  their  battle  against  the  Nan  Sung  in 
the  12th  year  of  Shih  Tsu  of  Yuan  (1275  A.D.)  to  burn  their  tents, 
etc.  The  "rocket",  "explosive  arrow"  and  "fire  arrow"  mentioned 
in  these  accounts  were  all  arrows  shot  out  to  set  fire  and  diff-~ 
from  the  real  rocket  which  is  propelled  by  the  force  generated  by 
igniting  explosives. 

Nevertheless,  one  should  realize  that  the  ancient  "rockets" 
would  lead  to  the  invention  of  the  true  rockets.  The  ancient 
"rocket"  was  thrust  out  by  the  elastic  force  in  the  bow  and  had 
igniting  power.  After  it  left  the  bowstring,  it  stopped  receving 
any  forward  push.  The  true  rocket,  however,  continues  to  be  acted 
upon  by  the  reaction  to  the  gas  jet  produced  by  combustion  during 
the  flight.  After  the  invention  of  explosives,  one  would  be  natur¬ 
ally  led  to  apply  explosives  or  fuses  to  the  arrow.  After  "rockets" 
had  been  used  for  a  long  time,  one  would  eventually  come  to  realize 
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that,  after  the  explosive  had  been  ifnli.^d.  It  rom/d.  i  no  would 
give  the  arrow  a  sustained  propulsion.  When  peep]*-  start. r-d  doing 
this  consciously,  the  rocket  was  born. 

II.  THE  INVENTION  OF  THE  ROCKET 

The  explosive  is  one  of  our  four  world  renowned  major  inven¬ 
tions.  It  is  a  necessary  factor  for  the  invention  of  the  rocket. 
With  the  equipment  of  explosives,  "rockets"  were  able  to  evolve 
into  the  real  rocket . 

Research  on  History  of  Chemistry  of  China  reveals  that  the 
"method  to  master  fire"  that  appeared  as  early  as  in  the  Tang 
Dynasty  employed  basically  the  three  ingredients  of  an  explosive, 
viz.  sodium  nitrate,  sulfur  and  carbon.  Three  recipes  [1]  for 
explosives  were  given  in  Wu  Ching  Tsung  Yao  compiled  by  Tsen  Kung 
Liang  of  Pei  Sung  (10*10  A.D.),  showing  that  explosives  had  been 
concocted  in  the  10th  Century  at  the  latest. 

The  rocket  was  probably  invented  some  time  between  the  10th 
and  11th  Centuries.  As  recorded,  in  the  third  year  of  Kai  Pao  of 
the  Sung  Dynasty  (970  A.D.),  "Defense  Department  Official  Feng  Chi 
Shen  et  al. suggested  the  method  of  rockets.  Order  was  given  to 
test  it  out  and  gifts  were  bestowed  of  them".  In  the  third  year 
of  Hsien  P’ing  (1000  A.D.),  "Ch'ang  T’ang  Fu  of  the  Shen  Wei  Shui 
troop  contributed  the  rockets,  fireballs  and  fire  caltrops  that  he 
made".  ("Ping  Chib  Ti  I  Pai  Wu  Shih",  vol.  197  of  Sung  Shih)  in 
the  above  records,  the  rockets  in  question  were  presented  to  the 
court  as  a  new  weapon.  Therefore,  these  could  not  be  the  "rockets" 
that  existed  in  the  Tang  Dynasty.  That  is  why  they  had  to  be 
tested  out  so  that  the  ruler  could  personally  witness  their  power 
and  effect.  If  these  were  of  the  old  type  of  "rockets",  then  they 
needed  not  to  be  tested.  The  T’ai  Tsu  and  T'ai  Tsung  brothers 
of  the  Sung  Dynasty  both  received  military  training  and  should  have 
been  familiar  with  that  type  of  "rocket".  Hence,  we  have  reason 
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to  bol  lew  that  ‘  bv  rockets  separately  offer,  d  by  Funr-  r'L  i  fhenr 
and  T' anr  Fu  were  the  precursors  of  the  i  ru'1  rocket ,  In  Vol .  12 


of  V/u  Chin-*  Tou.m"-  Yao  Ch '  1  f-*n  ('hi  is  recorded : 


To  set  off  the 


explosive: 


r.cur  five  ounces  of  explosives  into  thn  end  of 


the  arrow  to  which  birch  bark  had  been  added;  ignite  and  set  off". 
This  kind  of  "explosives  arrow"  was  probably  the  type  that  Fung 
and  T’ang  offered  and  probably  was  a  variation  of  the  true  rocket. 
The  following  premature  viewpoints  are  presented  to  invite  dis¬ 
cussion.  First,  it  was  mentioned  that  the  "explosives  arrow"  was 
to  be  "set  off".  The  words  "set-  off"  deserve  special  attention. 

They  can  only  apply  to  the  true  rocket.  For  the  "rocket"  of  the 
Tang  Dynasty,  the  appropriate  words  should  be  "shot  off"  instead 
of  "set  off".  Secondly,  based  on  the  wording  in  "ignite  and  set 
off",  one  can  see  that  the  dynamic  power  of  the  rocket  came  from 
ignition.  Obviously,  the  rocket  was  pushed  off  by  the  thrust  pro¬ 
duced  by  the  combustion  of  the  explosive.  Thirdly,  in  these 
accounts  no  mention  was  made  of  burning  the  opponent’s  tents.  From 
this,  it  is  apparent  that  the  rockets  were  not  the  "rockets"  used 
for  setting  fire.  Otherwise,  one  should  encounter  such  phrases  as 
to  "burn  down  the  enemy’s  tents".  Besides,  "pour  five  ounces  of 
explosives  into  the  end  of  the  arrow"  was  certainly  not  the  manner 
in  which  the  "bow-slung  fiery  pomegranate  arrow"  (Figure  1)  was 
loaded.  According  to  Huo  Lung  Chlng,  written  in  the  early  Ming 
Dynasty,  to  use  a  "bow-slung  fiery  pomegranate  arrow",  one  would 
"wrap  the  explosive  with  two  or  three  layers  of  tissue  paper, 
insert  the  shaft  of  the  arrow,  form  into  the  shape  of  a  pomegranate, 
wrap  tightly  with  burlap,  seal  with  melted  pine  resin,  ignite  the 
fuse  and  shoot  off  only  after  combustion  has  taken  place".  This 
is  different  from  "pouring"  In  the  explosive.  In  order  to  pour  in 
the  explosive,  one  would  need  a  cylindrical  object  of  sizable  open¬ 
ing  to  pour  into. 

The  three  points  mentioned  above  are  evidence  that  the  "explo- 
sives'arrow" ‘alluded  to  in  Wu  Ching  Tsung  Yao  Ch’ien  Chi  was  a  true 
rocket.  The  invention  of  the  rocket  probably  took  place  in  the 
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latter  part  of  t h o  10th  Century,  i.e.,  early  Pel  Hung.  From  the 
time  it  was  invented  to  the  time  it  was  included  in  books,  there 
ought  to  be  a  certain  time  interval.  The  1 1  m/*  lap  between  the 
time  Fung  Chi  Chen  offered  his  method  of  rockets  and  the  time  Wu 
Chine;  Tsung  Yao  was  completed  is  71  years,  and  that  between  the 
moment  T ' a n g  Fu  offered  his  rockets  and  the  completion  of  Wu  Ching 
Tsung  Yao  is  h 0  years.  It  is  logical  that'  only  after  the  value  of 
the  rocket  had  been  borne  out  by  actual  usage  in  battles  was  the 
rocket  included  in  books  on  military  strategy. 

Up  until  now,  most  people  believe  that  the  true  rocket  was 
invented  by  our  people  during  Chin  or  Nan  Sung.  Based  on  Chou 
Mi's  Wu  Lin  Chiou  Shih,  Liu  Hsien  Chou  estimated  the  time  to  be 
around  1250  A.D.  [2].  In  a  book  recently  published  in  Shanghai, 
it  is  said  that  the  rocket  was  "invented  in  the  latter  part  of  Nan 
Sung"  [3].  Somebody  even  said  vaguely  that  Wu  Pei  Chih  contains 
accounts  of  the  invention  of  rockets  by  our  people  before  the  Yuan 
and  Ming  Dynasties.  Europeans,  e.g.,  0.  G.  Sutton  and  Herbert  S. 
Zim,  believe  that  the  true  rocket  was  employed  by  the  Ching  people 
in  the  battle  being  Ching  and  Sung  in  1232  A.D.  [23.  We  regard 
these  statements  as  inaccurate  as  they  have  set  the  time  of  inven¬ 
tion  of  the  rocket  by  our  people  later  by  more  than  200  years. 
Japanese  Chi  Tien  Kuang  Pang,  in  a  composition  devoted  to  military 
tactics  of  the  Sung  and  Yuan  eras,  only  quoted  from  Sung  Huey  Yao 
the  account  of  T'ang  Fu  offering  the  rocket  and  did  not  explain 
what  kind  of  rocket  it  was  [53* 

In  the  past,  most  of  the  reference  material  used  by  our  scho¬ 
lars  came  from  Wu  Pei  Chih.  Although  its  accounts  of  the  rocket 
are  valuable,  they  were  given  too  late.  This  book  was  compiled  by 
Mao  Yuan  I  in  the  first  year  of  Tien  Hu  in  the  late  Ming  Dynasty 
(1612  A.D.).  Actually,  most  of  the  material  was  already  available 
in  the  early  Ming  Dynasty. 


Figure  1.  Bow-slung  fiery 
pomegranate  arrow 

(Chiao  Yu:  Huo  Lung  Chin^ 
vol.  2) 


Figure  2.  Rocket 

(Chiao  Yu:  Huo  Lung  Ching  vol,  2) 

l--assembled  rocket;  2 — dragon¬ 
shaped  stand;  3 — rocket 


The  aforementioned  Huo  Lung  Ching  was  written  in  the  10th 
year  of  Yung  Le  (1412  A.D.)  which  is  over  200  years  earlier  than 
Wu  Pei  Chih.  Much  of  the  material  and  many  figures  of  the  latter 
were  adopted  from  the  former.  The  author  of  Huo  Lung  Ching, 

Chiao  Yu,  was  a  specialist  of  weapons.  The  book  consists  of  three 
volumes.  Not  many  printed  copies  were  made  and  hence  are  bard  to 
find.  I  will  write  another  paper  specially  devoted  to  this  book. 

The  simplest  true  rocket  was  described  in  Huo  Lung  Ching  in 

words  and  with  figures  (Figure  2).  It  was  written  that  one  would 

"use  a  small  bamboo  stake  of  length  four  feet  and  two  inches,  an 

iron  arrowhead  4-1/2  inches  long  and  a  four  inch  long  iron  pendant 

to  be  attached  to  the  end.  Tie  a  short  cylinder  in  the  front  for 

making  fire.  When  ready  to  set  off,  place  in  a  dragon-shaped 

> 

stand  or  bamboo  container,  whichever  is  more  convenient".  Here 
again,  the  words  "set  off"  was  used  instead  of  "shoot  off".  That 
the  rockets  had  to  be  placed  in  a  dragon-shaped  stand"  or  "bamboo 
container"  supports  our  interpretation  of  the  accounts  in  Wu  Ching 
Tsung  Yao.  From  this  written  record,  one  can  see  that  there  already 
existed  a  uniform  specification  for  the  structure  of  the  rocket. 


42 


yr;- 


r 


FI ac Inr  the  rocket  in  a  container  or  on  a  .stand  helps  reduce 
stray  force  components  so  chat  the  rocket  can  bo  launched  along 
a  given  direction.  Moreover,  the  momentum  produced  by  the  com¬ 
bustion  of  the  explosive  will  not  be  dissipated,  thus  enabling 
the  rocket  to  fly  farther.  Rockets  of  such  design  could  not  have 
been  just  created,  but  must  have  bad  a  fairly  long  history  of 
development . 

III.  DEVELOPMENT  OP  THE  ROCKET 

The  Yuan  and  Ming  Dynasties  saw  a  great  deal  of  development 
of  the  rocket.  In  late  Yuan  and  early  Ming  at  the  latest,  a 
method  already  existed  for  launching  off  heavier  objects  by  means 
of  simultaneously  setting  off  many  rockets.  A  good  example  is 
found  in  the  "mysterious  fire  flying  duck"  recorded  in  Huo  Lung 
Ching  (Figure  3).  This  type  of  "flying  duck"  was  made  "with  bam¬ 
boo  strips,  with  size  of  a  duck  that  weighs  about  one  pound. 

Should  be  elongated  in  shape.  Use  paper  pulp  to  seal  and  shape. 
Fill  with  explosives  and  stuff  with  tissue  paper.  Attach  head 
and  tail.  Staple  paper  wings  at  the  sides  so  that  it  looks  like 
a  flying  duck.  Under  each  wing  place  two  big  lighting  sticks 
which  are  connected  to  fuses  brought  out  through  a  hole  in  the  back 
of  the  duck.  First  ignite  the  fuses.  After  the  duck  has  flown 
over  100  yards  and  is  about  to  drop  to  the  ground,  then  the  body 
starts  to  catch  fire  which  will  cause  a  big  fire  in  the  whole 
field".  The  method  was  thus  to  send  a  combustible  duck  to  the 
enemy's  side  by  means  of  many  rockets,  setting  fire.  As  far  as 
the  recording  goes,  the  same  figures  found  in  Huo  Lung  Ching  were 
also  found  in  another  book  by  Chiao  Yu,  Huo  Lung  Sben  Cb ' i  Chen  Fa 
(1377  A.D.)  which  dated  back  more  than  30  years  earlier. 

In  the  late  Yuan  and  early  Ming  Dynasties,  the  rocket  had 
another  development,  i.e.,  many  rockets  were  set  off  at  the  same 

I 

time.  According  to  Ming  Rhlh  Lu  Yung  Le ,  in  the  2nd  year  of  Chien 
Wen  (1400  A.D.),  the  so-called  "swarm  of  bees"  was  employed  in  a 
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Figure  4.  A  swarm  of  bees 

(Wu  Pei  Chih,  vol.  127) 

1 — disassembled  view;  2 — cover; 
3 — door;  4 — upper  arrow  rack; 

5 — lower  rack;  6 — 4  inch 
explosive  cylinder;  7 — 4  foot, 

2  inch  shaft 


Figure  5.  Fire  dragon  coming  out 
of  the  water. 

(Chiao  Yu:  Huo  Lung  Ching, 
vol.  3) 


battle.  This  was  a  powerful  weapon  which  could  "go  through  the 
bodies  of  both  man  and  horse".  A  "swarm  of  bees"  was  not  recorded 
in  Huo  Lung  Ching.  but  a  figure  of  it  is  found  in  Wu  Pei  Chih 
(Figure  4).  32  arrows  are  stored  in  a  wooden  container  and  can  be 

shot  out  to  a  distance  of  more  than  300  feet.  "After  the  main  fuse 
is  ignited,  all  the  arrows  are  set  off  simultaneously".  In  that 
book  were  recorded  many  other  similar  weapons,  some  with  several 
tens  of  arrows,  others  with  100  arrows,  etc. 


i t ■  vented  in  The  Yuan 


Sir;;: sin  t  wo-.;i  rockets  were  already 

and  Hin.r;  Dynasties.  A  repr^sentat  i  vn  in  Ui*.-  "  f  i  ro  dragon  eor.ing 

out  of  the  water"  (Figure  5)  given  in  vol.  3  of  Hue  Lunr  Ching. 

It  w an  recorued  thus:  "Take  a  five  foot  bamboo  stick,  remove  the 
nodes  and  scrape  it  thinner.  Use  a  wooden  carved  dragon  head  in 
the  front  and  carve  the  end  of  the  bamboo  stick  to  form  the  dra¬ 
gon's  tail.  Insert  several  arrows  in  the  dragon's  body  and  con¬ 

nect  all  the  fuses  to  an  outlet  in  the  dragon's  head.  When  fight¬ 
ing  on  water,  place  the  dragon  three  or  four  feet  above  the  water, 
ignite  and  it  will  fly  away  two  to  three  miles.  From  afar,  it 
looks  just  like  a  dragon  coming  out  of  the  water.  When  the  explo¬ 
sives  in  the  cylinder  arc-  almost  used  up,  the  arrows  fly  out  to 
burn  all  the  men  and  the  boat".  One  may  say  that  such  rockets  are 
fairly  advanced.  A  jet  propulsion  mechanism  was  used  in  the  spec¬ 
ially  fabricated  bamboo  cylinder  to  propel  the  entire  flight  of 
the  rocket.  From  Figure  5  one  can  see  that  four  rockets  were  used 
for  the  propulsion.  There  were  also  several  arrows  in  the  cylinder 
which  were  set  off  during  the  flight.  The  distance  of  flight 
reached  two  to  three  miles,  showing  that  there  was  a  fairly  large 
propelling  power. 

During  the  late  Ming  Dynasty,  the  rocket  became  even  more  ad¬ 
vanced — it  could  be  retrieved.  In  Wu  Pei  Chih  was  recorded  a  kind 
of  weapon  called  a  "flying  sand  cylinder"  (Figure  6).  "There  are 
various  ways  of  making  the  flying  sand  cylinder.  If  sand  from  the 
river  is  not  available,  one  can  crush  rocks,  sift  out  the  fine  dust 
with  loosely  woven  raw  silk  fabric  and  then  sift  out  the  sand.  To 
each  peck  of  sand  add  a  liter  of  explosives,  mix  over  heat  and 
set  aside.  The  mouth  of  the  firearm  is  made  of  thin  bamboo  slices. 
There  are  two  cylinders  containing  the  explosives.  These  are  to 
be  tied  together  in  such  a  manner  that  the  cylinder  in  the  front 
opens  to  the  back  while  the  cylinder  in  the  back  opens  to  the  front. 
Attach  to  the  cylinder  in  the  front  a  firecracker  seven  inches  long 
and  0.7  inches  in  diameter  which  is  wrapped  in  three  to  five  layers 


of  paper  and  glued  to  the  fire  start  "r .  *I’l i •. -  -rav'd  rand  arid 

explosive  mixture  is  placed  around  the  firecracker  and  tightly 
sealed.  Use  a  small  dart  at  the  ton.  If  the  rocket  is  to  be  used 
on  land,  omit  the  dart.  First  ignite  th-.  fir-.-  starter  that  will 
cause  the  rocket  to  move  forward,  using  a  big  bamboo  cylinder  as 
a  guiding  cylinder.  Aim  at  the  enemy’s  boat  and  the  dart  will  hit 
the  sail.  The  enemy  will  come  out  to  put  out  the  fire.  At  that 
moment,  the  first  cylinder  will  explode  and  the  sand  will  blind 
the  enemy.  The  second  cylinder  will  be  propelled  backward  to  our 
base  and  the  enemy  will  never  know  what,  actually  hit  them".  This 
is  also  a  primary  model  of  a  two-stage  rocket  but  is  even  better 
than  the  "mysterious  fire  flying  duck"  in  that  it  can  "be  propelled 
back  to  our  base" . 

From  the  above,  it  is  obvious  that  our  working  people  of 
ancient  times  had  many  inventions  in  the  field  of  the  rocket  and 
thus  had  many  glorious  accomplishments.  However,  in  the  latter 
part  of  the  l8th  Century,  especially  since  1840,  our  rocket  techno¬ 
logy  has  lagged  behind.  This  lag  is  entirely  due  to  the  oppression 
and  exploitation  inflicted  upon  us  by  foreign  capitalism  and  the 
reactionary  government  at  home.  In  the  semi-colonial,  semi-feudal 
old  China,  the  reactionary  forces  at  home  and  abroad  wildly  des¬ 
troyed  our  people's  economy  and  scientific  technology.  However, 
even  under  these  conditions,  our  working  people,  faced  with  the 
danger  of  the  destruction  of  our  nation  and  people,  still  kept  on 
working  on  various  weapons  including  the  rocket  in  the  hope  that 
they  might  use  these  to  fight  the  invaders.  In  189*1,  a  firearm 
maker  in  Chiang  Hsi  (anonymous)  "made  a  new  type  of  rocket,, the 
best  of  which  can  go  five  miles".  However,  the  government  at  that 
time  did  not  give  any  support  to  him.  The  test  flight  even  invoked 
displeasure  in  the  chancellor  [6]. 

A  new  turn,  however,  finally  took  place.  After  our  country  was 
liberated  in  19*19,  production  and  science  had  speedy  development 
under  the  wise  leadership  of  Chairman  Mao  and  the  Party.  Our  people 
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rockc-ts  and  satellites.  On  April  2 Jl , 
launch-.. a  our  first  manmade  satellite-  o 
of  1975,  we  further  accurately  landed 


too",  and  did  research  on 
1970,  wc  successfully 
f  the  Earth .  In  December 
it  on  the  Earth. 


At  present,  our  Party,  headed  by  Chairman  Hua,  has  summoned 
us  to  modernize  our  scientific  technolory-  As  a  result,  our  rocket 
and  space  technology  will  definitely  have  an  even  greater  develop¬ 
ment,  catch  up  with  and  surpass  world  standards,  and  make  yet 
greater  contributions  to  mankind. 


Figure  6.  Flying  sand  cylinder 
(Wu  Pei  Chih,  vol.  129) 

1 — bamboo  guiding  cylinder 
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